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ABSTRACT 

We compute the specific angular momentum distributions of a sample of low mass disk 
galaxies observed by Swaters. We compare these distributions to those of dark matter 
haloes obtained by Bullock et al. from high resolution TV-body simulations of struc- 
ture formation in a ACDM Universe. We find that although the disk mass fractions 
are significantly smaller than the Universal baryon fraction, the total specific angular 
momenta of the disks are in good agreement with those of dark matter haloes. This 
suggests that disks form out of only a small fraction of the available baryons, but yet 
manage to draw most of the available angular momentum. In addition we find that 
the angular momentum distributions of disks are clearly distinct from those of the 
dark matter; disks lack predominantly both low and high specific angular momentum. 
Understanding these findings in terms of a coherent picture for disk formation is chal- 
lenging. Cooling, feedback and stripping, which are the main mechanisms to explain 
the small disk mass fractions found, seem unable to simultaneously explain the disk's 
angular momentum distribution. In fact, it seems that the baryons that make up the 
disks must have been born out of angular momentum distributions that are clearly 
distinct from those of ACDM haloes. However, the dark and baryonic mass compo- 
nent experience the same tidal forces, and it is therefore expected that they should 
have similar angular momentum distributions. Therefore, understanding the angular 
momentum content of disk galaxies remains an important challenge for our picture of 
galaxy formation. 

Key words: galaxies: formation — galaxies: fundamental parameters — galaxies: 
kinematics and dynamics — galaxies: structure — dark matter. 



1 INTRODUCTION 



parameterize the angular momentum of dark matter haloes 
by the dimensionless spin parameter 



Disk galaxies are rotationally supported systems whose 
structure is governed by angular momentum. In the cur- 
rent paradigm of structure formation, galaxies form hier- 
archically by the assembly of dark matter haloes and the 
subsequent cooling of the baryonic mass component. In this 
standard picture, protogalaxies acquire angular momentum 
by tidal interactions with neighboring protogalaxies; the so- 
called cosmological torques (Hoyle 1953). This mechanism 
of spinning up the dark and baryonic matter has been stud- 
ied in great detail using numerical simulations (e.g., Barnes 
& Efstathiou 1987; Warren et al. 1992; Zurek, Quinn & 
Salmon 1988; Sugerman, Summers & Kamionkowski 2000) 
and found to be in good agreement with linear tidal-torque 
theory (e.g., Doroshkevich 1970; White 1984; Catelan & 
Theuns 1996). 




with A ~ 0.06 and a\ ~ 0.5 (e.g., Barnes & Efstathiou 1987; 
Ryden 1988; Cole & Lacey 1996; Warren et al. 1992). Start- 
ing with the seminal paper of Fall & Efstathiou (1980), it 
was soon realized that the size distribution of disk galaxies 
can be explained as originating from the spin parameter dis- 
tribution if the assumption is made that baryons conserve 
their specific angular momentum when cooling to form lu- 



where J, E, and M are the total angular momentum, energy 
and mass of the halo, and G is the gravitational constant. 
The distribution of A is well described by a log-normal dis- 
tribution, 




(2) 



Following Peebles (1969) it has become customary to 
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minous galaxies. This concept has spawn our current picture 
for the formation of disk galaxies, which has been addressed 
by numerous studies. Blumenthal et al. (1986) and Flores 
et al. (1993) investigated how adiabatic contraction of dark 
matter haloes impacts on the rotation curves of disk galax- 
ies. Kauffmann (1996) linked the formation of disk galaxies 
within this framework to the evolution of damped Lya ab- 
sorption systems. Dalcanton, Spergel & Summers (1997) and 
Mo, Mao & White (1998) investigated the structural prop- 
erties of disks, with emphasis on the variance induced by the 
A-distribution. Subsequent studies included recipes for bulge 
formation, gas viscosity, star formation and/or feedback and 
investigated more detailed properties of these model disk 
galaxies, such as the Tully-Fisher relation, the gas mass 
fractions, and the origin of the Hubble sequence (van den 
Bosch 1998, 2000, 2001; Jimenez et al. 1998; Natarajan 1999 
Heavens & Jimenez 1999; van den Bosch & Dalcanton 2000 
Firmani & Avila-Reese 2000; Avila-Reese & Firmani 2000 
Efstathiou 2000; Zhang & Wyse 2000; Buchalter, Jimenez 
& Kamionkowski 2001). 

The standard picture of disk formation that has 
emerged from these studies has been remarkably success- 
ful in explaining a wide variety of observational properties 
of disk galaxies. However, two important problems, both re- 
lated to the angular momentum of disk galaxies, have come 
to light. First of all, detailed hydro-dynamical simulations of 
this process of disk formation in a cold dark matter (CDM) 
Universe yield disks that are an order of magnitude too small 
(Navarro & Benz 1991; White & Navarro 1993; Steinmetz 
& Navarro 1999). This problem, known as the angular mo- 
mentum catastrophe, is a consequence of the hierarchical 
formation of galaxies which causes the baryons to lose a large 
fraction of their angular momentum to the dark matter. The 
second problem concerns the actual density distribution of 
disk galaxies. Under the assumption of detailed angular mo- 
mentum conservation, this density distribution is a direct 
reflection of the distribution of specific angular momentum 
of the protogalaxy. Although the distribution of total specific 
angular momentum of dark matter haloes is well established 
(i.e., equation little is known about the actual distribu- 
tion of specific angular momentum in each individual halo. 
Therefore, prompted by the observed surface brightness pro- 
files of disk galaxies, most models for the formation of disk 
galaxies have made the a priori assumption that the disks 
that form have exponential density distributions. Recently, 
however, Bullock et al. (2000) determined the angular mo- 
mentum distributions of individual dark matter haloes, and 
showed that disks that form out of such distributions are 
more centrally concentrated than an exponential. Van den 
Bosch (2001) has shown that when star formation, bulge for- 
mation, and feedback are taken into account, the resulting 
stellar density distributions are inconsistent with observa- 
tions, at least for the low surface brightness disk galaxies. 
In addition, van den Bosch showed that the truncation radii 
of the resulting disks are too small compared to observa- 
tions. 

From the above discussion it is clear that our under- 
standing of the formation of disk galaxies is directly related 
to the angular momentum distribution (hereafter AMD) of 
the baryonic mass component of protogalaxies. In this pa- 
per we use the observed density distributions and rotation 
curves of a sample of 14 dwarf galaxies to compute their 



AMDs. A comparison of these AMDs with those of dark 
matter haloes then provides important clues to the forma- 
tion process of disk galaxies. 

This paper is organized as follows. In Section |2| we 
discuss the angular momentum distributions expected on 
grounds of the standard paradigm of galaxy formation. In 
Section |^ we then present the angular momentum distribu- 
tions of a sample of 14 dwarf galaxies. The implications of 
these distributions for our picture of galaxy formation are 
discussed in Section ^, and we summarize our results in Sec- 
tion H 

Whenever the cosmological framework is important for 
our discussion we adhere to the currently popular ACDM 
model with fi = 0.3, tt A = 0.7, h = 0.7, and erg = 1.0. 
This model yields a reasonable fit to the current suite of 
cosmological constraints, including high redshift supernovae 
(Perlmutter et al. 1998; Riess et al. 1998; Garnavich et 
al. 1998), the cosmic microwave background radiation (e.g., 
de Bernardis et al. 2000), and the observed cluster abun- 
dances (Eke, Cole & Frenk 1996; Bahcall & Fan 1998). For 
the baryon density we adopt fibar = 0.019 hT 2 as suggested 
by the observations of primordial deuterium abundances by 
Tytler et al. (1999). With the cosmological parameters de- 
fined as above this implies a Universal baryonic mass frac- 
tion Of /bar = Obar/^0 = 0.13. 



2 THEORETICAL BACKGROUND 

Consider a virialized halo with mass M v i r consisting of dark 
and baryonic matter, and let p(j)dj indicate the fraction of 
mass with specific angular momentum between j and j + dj. 
The total specific angular momentum is given by 



Jtot — Jmax 



m(l)dl 



(3) 



Here jmax is the maximum specific angular momentum of 
P(j)> ' = j/imax, and m(j) is the normalized cumulative 
distribution 



m(j) 



P(j)dj. 



(4) 



The dimensionless spin parameter A can be related to jtot 
as 



A 



Jtot 



V^2 7*vir Vvi 



(5) 



Here r v ir is the virial radius, Kdr = \J GM v j r /r v ir, and Q is 
a geometrical factor. 

It is customary to parameterize the density distribution 
of virialized structures as 

^W=„(^)- 7 (l + ^) 7 - 3 (6) 

with r s a scale- length and p s = p v i r (r B ). For the density dis- 
tribution of the form ([]), which reduces to the NFW profile 
(Navarro, Frenk & White 1996, 1997) for 7 = 1.0, one finds 



g = G(l,C vir ) = VMCvir)/ _1 (Cvir), (7) 

where c v ir = r v i T /r s is the halo concentration parameter, 



f(x)= I dyy^il+yy- 3 , 



(8) 
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and 



h(x) 



Table 1. Properties of sample of late- type dwarf galaxies. 



dy/(»)y 1-7 (l + y) 



7-3 



(9) 



In the standard picture of disk formation it is assumed 
that baryonic and dark matter experience the same tidal 
torques, and thus end up with the same AMD p(j). As shown 
by Bullock et al. (2000, hereafter BOO) the AMD of dark 
matter haloes is well described by 



m(l) = 



fii 



1 + H-V 



(10) 



with I = j/jmax. For dark matter haloes in a ACDM Uni- 
verse with fi = 0.3, Q,a = 0.7, h = 0.7, and <j 8 = 1.0, BOO 
found that the distribution of /i is Gaussian in log(/i — 1) 
with a mean of —0.6 and a standard deviation of 0.4. Thus 
for 90 percent of the haloes 1.06 < fj. < 2.0 with a mean of 
(fi) = 1.25. For the probability distribution corresponding 
to ( p"o| ) one can write 

Cfi( P - 1) 



(Cs + M - 1) 
with s = j/jtot and 



Jtot 



= 1-At 



1 - (/x - l)ln 



/x-1 



(11) 



(12) 



If our current picture for the formation of disk galaxies is 
correct, the baryonic mass of disk galaxies thus should have 
an AMD of the form (p"l|). However, not necessarily all the 
baryons inside r v ir are currently part of the disk (stars plus 
cold gas). Part of the baryons may not have cooled by the 
present time, and/or a certain fraction of the baryons may 
have been expelled from the halo by either feedback pro- 
cesses (i.e., galactic winds) or by some stripping mechanism. 
We therefore define the fraction /disk = Mdisk/M v i r , where 
Mdisk is the total disk mass observed. If all available baryons 
have cooled to become part of the disk one expects that 
/disk = /bar- By comparing the AMD of the observed disk 
with that of equation (jll]), and by investigating how these 
AMDs correlate with /disk, important insights into the for- 
mation mechanism of dwarf galaxies can be obtained. 



3 THE ANGULAR MOMENTUM CONTENT 
OF DWARF GALAXIES 

3.1 Description of the data 

In order to compute the angular momentum distribution of 
the baryonic component of a disk galaxy one needs (i) the 
rotation velocity as function of radius, (ii) the distribution 
of the stellar mass, and (iii) the distribution of the gas. One 
thus requires a combination of both optical and HI observa- 
tions. Since most angular momentum is contained by mate- 
rial at large galactocentric radii it is essential that the obser- 
vations probe out to large enough radii. Recently, S waters 
(1999) obtained both optical (ii-band) and HI observations 
of a large sample of dwarf galaxies, down to fairly low HI 
surface brightness. These data are therefore ideally suited to 
compute angular momentum distributions. A subset of these 
data, with the highest quality rotation curves, was analyzed 
in detail by van den Bosch & Swaters (2001; hereafter BS01). 



UGC 


I) 


Mb 


f'o 


ft , 


* last 




731 


8.0 


— 16.63 


23.0 


1.65 


74 


57 


3371 


12.8 


— 17.74 


23.3 


3.09 


86 


19 


4325 


10.1 


— 18.10 


21.6 


1.63 


92 


11 


4499 


13.0 


—17.78 


21.5 


1.49 


74 


50 


6446 


12.0 


— 18.35 


21.4 


1.87 


80 


52 


7399 


8 1 


— 17.12 


20.7 


0.79 


109 


55 


7524 


3.5 


-18.14 


22.2 


2.58 


79 


16 


8490 


4.9 


-17.28 


20.5 


0.66 


78 


50 


9211 


12.6 


-16.21 


22.6 


1.32 


65 


14 


11707 


15.9 


-18.60 


23.1 


4.30 


100 


68 


11861 


25.1 


-20.79 


21.4 


6.06 


153 


50 


12060 


15.7 


-17.95 


21.6 


1.76 


74 


40 


12632 


6.9 


-17.14 


23.5 


2.57 


76 


46 


12732 


13.2 


-18.01 


22.4 


2.21 


98 


39 



Column (1) lists the UGC number of the galaxy. Columns (2) 
— (7) list the distance to the galaxy (in Mpc), absolute 7?-band 
magnitude, central J?-band surface brightness (in mag arcsec - 2 ), 
scale length of the stellar disk (in kpc), the observed rotation 
velocity Vi as t (in kms -1 ) at the last measured point, and the 
adopted inclination angle (in degrees), respectively. Magnitudes 
and central surface brightnesses have been corrected for inclina- 
tion and galactic extinction, but not for internal extinction. 

Taking account of the effects of beam smearing and adi- 
abatic contraction they determined the concentration and 
the virial velocities of the dark matter haloes, which they 
found to be consistent with predictions of ACDM models. 
In the following we concentrate on the 14 galaxies for which 
BS01 achieved a meaningful fit. These galaxies are listed in 
Table [j], together with some global properties taken from 
Swaters (1999). 



3.2 Mass modeling 

We now turn to computing the angular momentum distribu- 
tions of the 14 dwarf galaxies in our sample. For that purpose 
we first construct a mass model of each galaxy that best fits 
the observed rotation curve after the effects of beam smear- 
ing are taken into account. If the effects of beam smearing 
are small, one could in principle directly infer the angular 
momentum distribution from the observed V ro t(r) and the 
surface density distributions of the stars and gas. However, 
the advantage of the mass modeling is that typically HI is 
measured out to much larger radii than the maximum ra- 
dius out to which Vrot is determined (simply because a mea- 
surement of Ehi requires a lower S/N than Vrot). Since low 
surface brightness HI at large radii can still contain a signifi- 
cant fraction of the total mass, and has high specific angular 
momentum, it is essential that one takes this gas into ac- 
count. A proper mass model that fits the observed rotation 
curve makes predictions for the circular velocity at all radii, 
and can thus be used to compute a more complete angular 
momentum distribution. Furthermore, the mass modeling 
allows to correct for the effects of beam smearing. 

For our mass modeling we use the same method as em- 
ployed in BS01, which we briefly repeat here for complete- 
ness. We assume that there are three mass components in 
each galaxy: an infinitesimally thin gas disk, a thick stellar 
disk, and a spherical dark matter halo. 

We make the assumption that the gas is distributed 
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Figure 1. The normalized cumulative angular momentum distribution, m(j), as function of j/jmax- Each panel plots, with thick lines, 
three distributions for a particular sample galaxy (as indicated in the panel) that correspond to three different values of the stellar 
mass-to-light ratio: T^j = (dotted lines), T^j = 1.0 (M/L)q (solid lines), and T^j = 2.0 (M/L)q (dashed lines). The three thin lines 
correspond to AMDs of the form ( Jioj ) with /i = 1.06 (upper curves), fj, = 1.25 (middle curves) and fi = 2.0 (lower curves) and are plotted 
for comparison. These values of fj, correspond to the mean and the 90 percent range of the distribution in /i found by Bullock et al. (2000) 
for ACDM haloes. Note that the majority of the disks has an AMD that is clearly distinct from that of the dark matter. 



axisymmetrically in an infinitesimally thin disk, and model 
the HI density distribution as: 

EhiCR) = £ (f-)^-* 7 * 1 +/Eoe-««-«^ 2 . (13) 

The first term represents an exponential disk with scale 
length R\ and with a central hole, the extent of which de- 
pends on p. The second term corresponds to a Gaussian ring 
with radius R2 and a FWHM oc a. The flux ratio between 



these two components is set by /. The form of equation ( |13| ) 
has no particular physical motivation, but is an appropriate 
fitting function, which yields excellent fits to the observed 
HI surface brightness distributions (see BS01)Q When com- 

* In the case of UGC 7524 the fitting function of equation 
can not satisfactorily describe the data and linear interpolation 
between the data points is used to describe £ hi (R). 
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Figure 2. Histograms (hatched) ofthe distribution ofAdisk f or t he 14 dwarf galaxies in our sample. These values have been determined 
from the AMDs plotted in Figure M using equations (31) and (|l7|). Results are plotted for three values of the stellar mass-to-light ratio, 
as indicated in each panel. The thick solid lines plots the probability distribution of equation (0) with A = 0.06 and o\ = 0.5, and is to 
represent the spin parameter distribution of dark matter haloes. Especially for Yjj = 1.0 (M/L)q the two distributions agree remarkably 
well, suggesting that the disks have the same distribution of total specific angular momentum as the dark matter haloes. 



Table 2. Results for T R = 1.0 (M/L) and 7 = 1.0. 



UGC 


Cvir 


v viI 


/"disk 


/gas 


"^disk 


Jtot 


J max 


731 


17.6 


48.6 


0.024 


0.801 


0.061 


308 


775 


3371 


10.6 


65.5 


0.018 


0.715 


0.056 


569 


1618 


4325 


33.0 


48.6 


0.037 


0.541 


0.074 


328 


971 


4499 


2.4 


126.3 


0.003 


0.672 


0.007 


330 


1195 


6446 


9.1 


56.2 


0.041 


0.570 


0.052 


397 


1325 


7399 


19.9 


65.8 


0.012 


0.691 


0.044 


396 


1692 


7524 


6.4 


78.8 


0.012 


0.500 


0.025 


393 


1025 


8490 


17.5 


53.2 


0.026 


0.769 


0.062 


378 


1106 


9211 


19.2 


41.2 


0.055 


0.865 


0.107 


381 


1058 


11707 


14.6 


62.2 


0.062 


0.770 


0.103 


886 


2046 


11861 


16.4 


93.1 


0.068 


0.405 


0.099 


1861 


4820 


12060 


31.1 


42.8 


0.102 


0.710 


0.168 


582 


1477 


12632 


16.5 


47.8 


0.033 


0.760 


0.078 


387 


976 


12732 


9.0 


68.9 


0.040 


0.869 


0.081 


926 


2267 



Column (1) lists the UGC number of the galaxy. Columns (2) 
and (3) list c vlr and V v i T (in kms - 1 ) ofthe best-fit mass model, 
respectively. Columns (4) and (5) list the disk mass fraction /disk 
and the gas mass fraction / ga s, respectively. Column (6) lists the 
baryonic spin parameter Adi s k> and columns (7) and (8), list the 
total and maximum specific angular momentum of the baryons 
(both in kpckms" 1 ). 

puting the circular velocities of the gas, we multiply Ehi by 
a factor 1.3 to correct for the contribution of Helium. 
For the stellar disk we assume a thick exponential 

p*(R, z) = pi exp(-R/R d ) sech 2 (z/zo) (14) 

where Rd is the scale length of the disk in the i?-band. 
Throughout we set zo = Rd/G. The exact value of this ratio, 
however, does not significantly influence the results. None of 
the galaxies in our sample has a significant bulge component. 

For the DM component we consider a density distribu- 
tion of the form ([]). Unless stated otherwise we focus on 
haloes with 7 = 1.0 (i.e., NFW profiles). Rather than pa- 
rameterizing the DM halo by (C200, V200), as in BS01, we use 
(cvir, Kir). Here 



Cvlr = I^ = C200 ^_|j , (is) 

with A v i r the virial density, defined as the average density 
inside r v i r expressed in terms of the critical density for clo- 
sure. For the ACDM cosmology used here A v i r — 101 (e.g., 
Bryan & Norman 1998). 

When fitting the rotation curves we take adiabatic con- 
traction (Barnes & White 1984; Blumenthal et al. 1986; Flo- 
res et al. 1993) and beam smearing into account (see BS01 
for details). The best fit values for c v i r and V v ii are listed 
in Table H together with the inferred values for /disk and 
/ gas . Note that the values of c v i r and Vvir differ somewhat 
from the values of c and V200 quoted in BS01 owing to the 
different definitions. 

3.3 Angular momentum distributions 

After finding the mass model that best fits the observed 
ShiM an d Vrot(r) we determine V c (r), L(r), and Mhi(t*) on 
a linear radial grid between r — and r = r max . Here V c (r) 
is the total circular velocity at r of the best fit model, L(r) 
is the total .R-band luminosity inside r, as computed from 
the thick exponential, Mhi(?") is the total HI mass inside r 
computed from the best fit HI surface density distribution, 
and r max is the maximum radius out to which HI is detected. 
We define the total disk mass inside radius r as 

A/disk(r) = T R ■ L(r) + 1.3 ■ M m {r). (16) 

We thus assume that the stellar mass-to-light ratio is con- 
stant with radius, and that the total gas mass is simply 1.3 
times the HI mass to take account of Helium. Any contribu- 
tion from molecular and/or ionized gas is therefore consid- 
ered negligible. The distribution of specific angular momen- 
tum is simply given by m(j) = M disk (r)/M disk (r ma , x ) with 
j = rV c (r). 

In Figure hi we plot m(j) as function of j'/jmax for all the 
galaxies in our sample. Results are shown for three different 
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Figure 3. The upper left panel plots c v i r versus V vlI for the mass models with Tr = 1.0 (M/L)q that best fit the observed rotation 
curves. The solid and dashed lines correspond to the median and 2<r limits of the distribution of c v ; r for the ACDM cosmology adopted 
here (computed using the model of Bullock et al. 2001). Although the dwarf galaxies reveal a somewhat steeper decline of c v ; r with 
increasing V v i T , reflecting the degeneracies in the rotation curve fitting, the overall distribution of halo concentrations is in reasonable 
agreement with the ACDM cosmology (see BS01 for a more detailed discussion). The upper right panel reveals a narrow correlation 
between Ab ar (determined from the AMDs) and /b ar (inferred from the best fit mass models). The implications of this correlation for 
the formation of disk galaxies are discussed in detail in the text. The lower left panel plots log(j'max) of the disk versus log(j'tot). The 
strong correlation found implies a narrow distribution of f = jtot/jmaxi which is shown in the lower right panel (shaded histogram). The 
thick solid curve corresponds to the distribution of £ of dark matter haloes in a ACDM cosmology, and is broader and offset to lower 
values of f . This reiterates that the AMDs of the disks in our sample are different from the AMDs of dark matter haloes found by BOO 
(see also Figures El and 0) . 



values of the i?-band mass-to-light ratio: = (dotted 
lines), T R = l.O(M/L) (solid lines), and T R = 2.O(M/L) 
(dashed lines). For comparison, we also plot the specific an- 
gular momentum distribution of equation (Jl^) for three val- 
ues of /i: 1.06, 1.25, and 2.0 (thin solid lines). These dis- 
tributions outline the 90 percent interval of the AMDs of 



ACDM haloes (see Section |2|). As is immediately apparent, 
the AMDs of the disks are clearly different from those of 
ACDM haloes. 

From the AMDs thus obtained we compute the total 
specific angular momentum of the disk, j to t, (using equa- 
tion which we parameterize by 
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Table 3. Correlation matrix. 





jtot 


J max 


/disk 


/gas 


Kir 


^ vir 


Adisk 


0.398 


0.231 


0.886 


0.407 


-0.543 


0.464 


jtot 




0.886 


0.538 


-0.055 


0.332 


-0.262 


J max 






0.363 


-0.024 


0.516 


-0.270 


/disk 








0.134 


-0.411 


0.345 


/gas 










-0.446 


-0.099 


Vvir 












-0.701 



Matrix of Spearman Rank's correlation coefficients, r s , for various 
parameters of the dwarf galaxies in our sample. A value of r s = 
+1.0 (—1.0) indicates a perfect correlation (anti-correlation), 
whereas values of r s near zero indicate that the parameters are 
uncorrelated. The only three significant correlations are those be- 
tween (A diak ,/ diak ), (jtotjmax), and (V v ir,Cvir)- For all other sets 
of parameters the values of r s obtained are consistent with no 
correlation. 



Adisk — 



Jtot 



■5(1-0, Cvir 



(17) 



y2 7*vir Vvi 

The values of jtot and Adisk are listed in Table ^| In Figure ^ 
we plot histograms of the normalized probability distribu- 
tions of Adisk for all three values of T_r. For comparison the 
probability distribution of equation (Q) with A = 0.06 and 
a\ = 0.5 is plotted as thick solid lines. As can be seen, the 
spin parameter distributions of the disks are fairly similar to 
that of the dark matter haloes, especially for realistic mass- 
to-light ratios in the range 1.0 (M/L)© <, T R <, 2.0 (M/L)©. 
This implies that, on average, the total specific angular mo- 
mentum of the disks is similar to that of the dark matter. 
Yet, as we have seen from Figure jj], their distributions of 
specific angular momentum are clearly distinct from that 
of the dark matter haloes. Furthermore, the baryonic mass 
fractions inferred from the best fit models to the observed 
rotation curves are significantly smaller than the Universal 
value /bar- All this holds important clues to the formation 
of (dwarf) galaxies, which we discuss in the next section. 



4 THE FORMATION OF DWARF GALAXIES 

In what follows we restrict ourselves to the AMDs for T_r = 
1.0 (M/L)©, which is a realistic value for the stellar mass- 
to-light ratio given the typical colors of the dwarfs in our 
sample (B S01 ; Bell & de Jong 2001). As we address shortly 
in Section 4.3 below, none of our results are sensitive to this 



particular value of T_r. 



4.1 Clues from parameter correlations 

In order to assess whether any significant correlations exist 
among the various parameters listed in Table ^| we compute 
the Spearman rank coefficients r s . The resulting correlation 
matrix is listed in Table ^. The only significant correlations 
are those between Vvir and c v ir (the probability of obtaining 
r s = —0.701 under the null hypothesis that no correlation 
exists is p s — 5.2 x 10 -3 ), and between Adisk and /disk, and 
/max and jtot (both with r s — 0.886 andp s = 2.5 x 10~ 5 ). We 
now address the meaning of each of these three correlations, 
which are plotted in Figure Q 

The decrease of c v ir with increasing Vvi r is, to first order, 
consistent with what one expects for a ACDM cosmology. 
This is indicated by the solid and dashed lines in the upper 



right panel of Figure which correspond to the mean and 
the 2a limits of the distribution of halo concentrations as 
predicted by the Bullock et al. (2001) model for the ACDM 
cosmology adopted here. However, given the amount of scat- 
ter expected, it is surprising that we find such a strong anti- 
correlation; i.e., the slope of the relation between Vvir and 
c v ir of the dwarf galaxies in our sample seems too steep. 
This most likely reflects degeneracies in the rotation curve 
fitting: a small under-(over)estimate of Vvir can result in a 
relatively large over-(under)estimate of c v i r (see BS01 for a 
more detailed discussion). 

The strong correlation between jtot and j ma x signals a 
small dispersion in the distribution of f. This is shown in 
the lower right panel of Figure [| The hatched histogram 
corresponds to the distribution in £ for the 14 dwarf galax- 
ies in our sample. We can compare this to the (^-distribution 
expected for dark matter haloes by converting the distribu- 
tion in n found by BOO using the correspondence between fx 
and £ (equation fl2|). The resulting probability distribution 
p(£) is plotted as a thick solid line. Comparing p(£) of the 
dark matter haloes with that of the dwarf galaxies one notes 
two important differences. First of all, the latter is signifi- 
cantly narrower than the former (which explains the strong 
correlation between j to t and j ma x)- Secondly, the mean ( of 
the dwarf galaxies is larger than that of the dark matter 
haloes. Low values of £ imply AMDs with a long tail to high 
specific angular momentum, and such AMDs are apparently 
underrepresented in dis k gal axies compared to dark matter 
haloes (but see Section 4.3.2 below). 

The correlation between Adisk and /disk, in the sense 
that dwarf galaxies with a higher inferred disk mass frac- 
tion have a higher specific angular momentum, is remark- 
ably strong. Based on the same set of data, but using a 
less detailed computation of Adisk based on the stellar com- 
ponent only, a similarly strong correlation between Adisk 
and /disk was recently found by Burkert (2000b). At first 
sight such correlation is consistent with an inside-out forma- 
tion scenario: the baryons that make up stars and cold gas 
have cooled from the inside out to form the galaxies. Since 
most angular momentum is contained in the outermost mass 
shells, which cool latest, one thus expects a correlation as 
seen. However, the problem with this interpretation is that if 
only a small fraction of the baryons have cooled, one would 
expect the total specific angular momentum of the disks to 
be significantly smaller than that of the dark matter. Yet, 
the distribution of A d i s k is comparable to that of dark mat- 
ter haloes, whereas /disk is significantly smaller than /bar- 
Apparently, galactic disks in low-mass systems form only 
out of a small fraction of the total baryonic mass, but yet 
manage to draw most of the available angular momentum. 
This puzzling aspect has recently also been noted by Navarro 
& Steinmetz (2000), using simple scaling relations of disk 
galaxies. 



4.2 Challenges for the standard model of disk 
formation 

In the standard picture, disks form out of baryons that cool 
and conserve their specific angular momentum. Since both 
baryons and dark matter experience the same tidal forces, it 
is generally assumed that both have the same angular mo- 
mentum distribution. Therefore, initially the AMD of the 
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s= j/jtot 

Figure 4. The shadeclareas indicate the AMDs p(s) of the 14 disk galaxies in our sample, normalized to /disk//bar- For comparison we 
plot p(s) of equation (nil) with [i = 1.25 (normalized to unity), and which represents the median of the AMDs of ACDM haloes. Under 
the standard assumption that baryons conserve their specific angular momentum the difference between the two distributions reflects 
the AMD of the baryonic matter that is not incorporated in the disk. Note that it is preferentially the baryonic matter with both the 
highest and the lowest angular momentum that is absent in the disks. 



baryonic component, with mass /bar-Wvir, can be parame- 
terized by equation (p"l]). A comparison of this distribution 
with the probability distributions p(s) of real galaxies can 
be used to test this standard picture and/or provide use- 
ful insights regarding the details of galaxy formation. Such 
comparison is only valid if the total baryonic mass and the 
disk material have the same j'tot (after all s = j/jtot). The 
fact that the distribution of Adisk for our sample of dwarf 
galaxies is consistent with p(X) of dark matter haloes (cf. 



Figure |2J) shows that this is a valid assumption to make, at 
least in a statistical sense. 

In Figure ^ we plot the distributions p(s), normalized 
to /disk//bar, for all galaxies in our sample (hatched areas). 
In addition we plot p(s) of equation ( |TT| ) with fi = 1.25 
(the median value found by BOO) and normalized to unity 
(solid lines) . This distribution is to represent the AMD of the 
total baryonic mass. Then, under the standard assumption 
that baryons conserve their specific angular momentum, the 
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difference between the two distributions describes the AMD 
of the baryonic material that did not make it into the disk. 

The first characteristic to note from Figure ^ is that the 
area under the hatched curves is much smaller than that un- 
der the solid lines, reflecting the fact that /disk is significantly 
smaller than the Universal baryon fraction /b ar . Secondly, in 
all cases p(s) of the total baryonic matter extends to much 
higher specific angular momenta. This indicates that its ra- 
tio Jmax/Jtot is larger than that of the disk material (cf. lower 
right panel of Figure ^) , and implies that none of the bary- 
onic material with the highest specific angular momentum 
has made it into the disk component. Finally, the difference 
between p(s) of the total baryonic mass and that of the disk 
mass is largest for small s. 

We can thus conclude that if the assumptions made, and 
which reflect our standard picture for the formation of disk 
galaxies, are correct, the combined effects of cooling, feed- 
back, and stripping have to be such that: (i) only a small 
fraction of the available baryons end up in the disk compo- 
nent, (ii) the total specific angular momentum material of 
the disk is comparable to that of the total baryonic mass, 
(iii) none of the highest specific angular momentum mate- 
rial makes it into the disk, and (iv) there is a tendency to 
preferentially prevent low angular momentum material from 
ending up in the disk. 

4.3 Uncertainties in the angular momentum 
distributions 

Before drawing any conclusions regarding the formation of 
(dwarf) disk galaxies, it is worthwhile to examine some of 
the uncertainties in the disk AMDs derived above. 



4.3.1 The Mass of the Disk 

One potential worry is that we have missed a significant frac- 
tion of the actual disk mass. In the discussion above we have 
focussed on results for Tr = 1.0 (M/L)©, which, for a Scalo 
(1986) IMF, yields colors in good agreement with the data 
(see discussion in BS01). However, internal dust extinction 
could have reddened the galaxies and/or a large fraction of 
brown dwarfs could be present. Furthermore, a central con- 
centration of molecular gas could be present, which would 
add low angular momentum material. These uncertainties in 
the disk's mass distribution can be modeled by varying Tr, 
either as a constant (mimicking for example the presence of 
brown dwarfs), or as function of radius (mimicking for exam- 
ple the presence of molecular gas). The freedom in Tr(R) 
is limited by the requirement that the mass model has to 
yield a reasonable fit to the observed rotation curve. We 
have performed extensive tests with varying Tr(R). Within 
the restrictions imposed by the rotation curves we find that, 
although it is possible to bring p(s) of the disk in some- 
what better agreement with that of the total baryonic mass, 
the four characteristics listed above remain. This is illus- 
trated in the upper panels of Figure |^ where we plot p(s) 
of UGC 12632 for three values of Tr (taken to be constant 
with radius). Even for Tr = 8.0 (M/L)©, which is the max- 
imum mass-to-light ratio allowed by the observed rotation 
curve, there is still a pronounced deficit of low-angular mo- 
mentum material in the disk compared to the total baryonic 
mass component. 



4.3.2 The Size of the Disk 

The maximum specific angular momentum of the disk, j max , 
is directly proportional to the size of the disk r max (rotation 
curv es a re generally flat at large radii). As stated in Sec- 
tion 3.3, we assume that r max is equal to the radius out to 
which HI is detected. Typically, these radii coincide with 
HI column densities of the order of ~ 10 19 cm~ 2 . This is 
close to the column density below which one expects the 
gas to be ionized by the cosmic background flux of ionizing 
photons (e.g., Sunyaev 1969; Maloney 1993). Therefore, it 
might well be that the actual gas disk extends significantly 
beyond the r max adopted here. If we make the assumption 
that this gas follows the same surface density profile as the 
HI we can compute the total disk mass missed by integrat- 
ing equation ( |13[ ) out to infinity, and comparing that to the 
total gas mass inside r max . We find that for all galaxies in 
our sample the gas at radii beyond r max does not contribute 
more than 0.5 percent of the total disk mass. Thus, whereas 
the uncertain extent of the disk can imply values of j max 
that are significantly underestimated, this does not influ- 
ence our estimates of jtot- For instance, if the actual size 
of the disks is about 1.4 times r max , the mean of £ would 
be in much better agreement with the expected mean for 
dark matter haloes (lower right panel of Figure^]). Also, the 
thick curves in Figure |l| would shift to the right by about the 
same factor, bringing them in somewhat better agreement 
with the AMDs of BOO. However, Figures ^ and ^ would 
remain virtually unchanged, and except for aspect (iii), the 
main problems outlined above thus remain. 



4-3.3 Asymmetric Drift 

Another problem with the AMDs derived in Section ^ is re- 
lated to the fact that we have assumed that the stars move 
on purely circular orbits. In reality, however, stars have a 
non-zero asymmetric drift, and we have thus overestimated 
the angular momentum of the stellar component. In order 
to assess the importance of asymmetric drift we perform 
the following test. We define a free parameter rj that de- 
scribes the (constant) ratio between the true rotation ve- 
locity of the stars and the local circular velocity]^]. In the 
lower three panels of Figure ^ we plot p(s) of UGC 12632 
for Tr = 1.0 (M/L)q and three different values of rj: 0.2, 
0.6, and 1.0 (i.e., no asymmetric drift). Only for rj £ 0.3 
do we find a discernibly different AMD. Typical asymmetric 
drifts correspond to values of rj much closer to unity, and we 
therefore conclude that our assumption that r\ — 1.0 does 
not significantly influence our results. This can also be un- 
derstood directly from the high gas mass fractions found for 
the dwarf galaxies in our sample. With 0.4 / gas 0.9 
the disk mass, and thus the AMD, of our dwarf galaxies is 
generally dominated by the gas. A small shift of the angu- 
lar momentum distribution of the stars relative to the gas 
therefore does not significantly influence the total AMD. 



t With this definition the asymmetric drift velocity is given by 
v a (r) = (l-v)V c (r). 
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Figure 5. The influence of mass-to-light ratio and asymmetric drift on the AMD of UGC 12632. The upper panels plot p(s) of 
UGC 12632 for three different values of the stellar mass-to-light ratio (as indicated in each panel). Since the stellar disk is more centrally 
concentrated than the HI component, increasing preferentially increases the disk mass fraction at low angular momentum. However, 
even for = 8.0 (M/L)q, the maximum mass-to-light ratio allowed by the rotation curve, a clear deficit at low s remains. The lower 
panels plot p(s) but now for three different values of n (as indicated in each panel), defined as the ratio of the rotation velocity of the stars 
to the circular velocity (r\ is thus a measure of the asymmetric drift). Only for the extreme asymmetric drift corrections with n = 0.2, 
does the AMD of UGC 12632 change its characteristic shape. The solid lines are as in Figure W 



4-3.4 Adiabatic Contraction 



4.4 Cooling, stripping, feedback & viscosity 



In the mass modelling used to fit the observed rotation 
curves we have made the assumption that the density dis- 
tribution of the dark matter halo is modified by adiabatic 
contraction. Since the assumptions underlying the method 
used (i.e., sphericity, adiabatic invariance) are not necessar- 
ily accurate it is wortwhile to examine what effect this adi- 
abatic contraction has on our results. To test this we have 
reanalyzed all data without any correction for adiabatic con- 
traction. The effects, however, are negligible. Typically one 
finds values of c v ir that are slightly larger (the dark matter 
halo now needs to be more centrally concentrated in order 
to fit the observed rotation curve), but the actual AMDs 
remain virtually unchanged. One parameter that is directly 
influenced by the change in c v i r is the spin parameter A 
(equation jH)), but even here the differences are small (less 
than 5 percent). 

Thus, despite considerable freedom in the stellar mass- 
to-light ratios, the sizes of the disks, and the asymmetric 
drift corrections, we conclude that the characteristics of the 
AMDs of disks outlined above are robust. As for the param- 
eters listed in Table ^: / gas , Adisk and j to t are very robust, 
Jmax may be significantly underestimated, and the uncer- 
tainties on Cvir, Vvir, and / d i s k are strongly correlated and 
discussed in more detail in BSOT Typically, lowering V v i T 
implies larger c v i r and /disk- This, through equation (g) re- 
sults in a slightly larger Adisk- As long as the quality of the 
rotation curve fit remains similar, the total specific angu- 
lar momentum, however, remains virtually unchanged. This 
implies that of the four main characteristics listed in Sec- 
tion 4.2 above, only (iii) is somewhat questionable as it is 



If cooling in small mass haloes is very inefficient, it might 
explain aspects (i) and (iii). After all, cooling is an inside out 
process, and the highest angular momentum material resides 
in the outs kirt s of the haloes. However, as already mentioned 
in Section 4.1, this picture is inconsistent with (ii), as one 



directly related to the assumed size of the disk. 



would expect Adisk to be much smaller than the spin param- 
eter of the dark matter haloes, and it does not explain aspect 
(iv). Furthermore, cooling is supposed to be extremely effi- 
cient in dense low mass haloes, and one typically expects the 
vast majority of the baryons in dwarf galaxies to have cooled 
by the present time. Tidal stripping has virtually the same 
effects as inefficient cooling; stripping is most efficient near 
the outskirts of the haloes, and thus might explain both (i) 
and (iii). However, as with cooling, stripping is inconsistent 
with both (ii) and (iv). 

The problem with feedback as the dominant cause for 
the AMDs observed is that one would expect a relatively 
strong correlation between /disk and Vvir. After all, the effi- 
ciency for blowing a galactic wind should be inversely pro- 
portional to the square of the galaxy's escape velocity, which 
in turn is proportional to Vvir- However, the correlation 
between /disk and Vvir has a Spearman rank coefficient of 
r s — —0.411 and thus seems to indicate an anti-correlation 
(although with a probability of p 3 = 0.14 this is not signif- 
icant). Furthermore, as shown by van den Bosch (2001), a 
simple model for supernovae induced galactic winds does not 
reveal any tendency to preferentially expel the low angular 
momentum material as required. 

One possibility is that one of the standard assumptions 
is wrong and that baryons do not conserve their specific an- 
gular momentum. However, given that disks have the same 



© 2000 RAS, MNRAS 000, 



The Angular Momentum Content of Dwarf Galaxies 11 



distribution of j to t as the dark matter haloes implies that 
one can only consider a scenario in which the specific an- 
gular momentum of the baryons is redistributed. However, 
the natural mechanism for angular momentum redistribu- 
tion, viscosity, transports low angular momentum inwards 
and high angular momentum material outwards. Therefore, 
any viscosity present in the disk will only aggravate the dis- 
agreement at low and high s between the AMDs of the disk 
and the total baryonic matter. 

4.5 The nature of dark matter 

One interesting interpretation of the large differences be- 
tween the AMDs of disks and cold dark matter haloes 
found might be that dark matter is not cold. Recently, 
several studies have focussed on warm and self-interacting 
dark matter (hereafter WDM and SIDM, respectively) sce- 
narios, and have shown that this results in dark matter 
haloes with constant density cores and with less substruc- 
ture (e.g., Spergel & Steinhardt 2000; Burkert 2000a; Colin, 
Avila-Reese & Valenzuela 2000; Bode, Ostriker & Turok 
2000; Dave et al. 2001; Kochanek & White 2000; Moore 
et al. 2000; Yoshida et al. 2000a,b). Not only might this 
solve the angular momentum catastrophe mentioned in Sec- 
tion [j] (Sommer-Larsen & Dolgov 2001), it might also allevi- 
ate the potential problems with the rotation curves of dwarf 
and low surface brightness galaxies (Flores & Primack 1994; 
Moore 1994; Burkert 1995; McGaugh & de Blok 1998; van 
den Bosch et al. 2000; van den Bosch & Swaters 2001) and 
with the abundance of satellite galaxies (Kauffmann, White 
& Guiderdoni 1993; Moore et al. 1999; Klypin et al. 1999). 
It would be worthwhile to examine to what extent these dif- 
ferences in halo structure result in AMDs that are different 
from those of CDM haloes. For instance, Moore et al. (2000) 
have pointed out that the central regions of SIDM haloes are 
rotationally supported, implying AMDs with less low angu- 
lar momentum material compared to CDM haloes.. Detailed 
high resolution iV-body simulations of structure formation 
in WDM and SIDM cosmologies are required to test whether 
a modification of the nature of dark matter brings the AMDs 
of disks and dark matter haloes in better agreement. 



4.6 Decoupling the baryons from the dark matter 

Probably the most likely conclusion from our results is that 
the standard assumption that baryons and dark matter have 
the same AMD is incorrect. Since dark and baryonic mat- 
ter should in principle experience the same cosmological 
torques, we thus need a mechanism to somehow decouple 
the baryonic mass component from that of the dark matter. 
One such mechanism, suggested by Hogan (1979), is Comp- 
ton drag on the background radiation. This causes the gas to 
be frozen into the comoving frame such that it experiences 
no tidal torques. However, this mechanism is only efficient at 
very high redshifts, and will cause the baryons to have less 
specific angular momentum than the dark matter. This is 
inconsistent with what is required, since, as shown by Burk- 
ert (2000b), the strong correlation between /disk and Adisk 
found requires a mechanism that can actually spin-up the 
baryons relative to the dark matter. 

A more plausible mechanism might be to resort to feed- 



back from early structure formation which stirs up the bary- 
onic component resulting in a density distribution that is 
quite distinct from that of the dark matter. Subsequent 
torques will then be different for the two mass components, 
causing a decoupling of their angular momentum distri- 
butions. Furthermore, it is important to realize that be- 
cause of the dissipative nature of the baryonic matter, an- 
gular momentum redistribution in merging systems will be 
quite different for baryonic and dark matter. High resolu- 
tion hydro-dynamical simulations should be useful to inves- 
tigate whether the AMDs of baryons and dark matter are 
similar (as usually assumed), or whether the processes men- 
tioned above cause a decoupling that bring the AMDs of the 
baryonic component in better agreement with those of disk 
galaxies. 



5 SUMMARY 

In the standard picture of structure formation the angu- 
lar momentum of protogalaxies originates from cosmologi- 
cal torques. Since dark and baryonic matter experience the 
same tidal forces, it is expected that both mass components 
end up with the same distribution of specific angular mo- 
mentum. In a recent paper Bullock et al. (2000) determined 
the AMDs of dark matter haloes in a ACDM cosmology. If 
our picture of disk formation is correct, we thus expect disks 
to be born out of similar distributions. A comparison of the 
AMDs of disks with those of dark matter haloes thus yields 
important insights into the formation mechanism of galax- 
ies, and may be used to test our standard picture of disk 
formation. 

In this paper we have computed the AMDs of a sample 
of 14 dwarf galaxies with observed HI rotation curves. These 
rotation curves have been analyzed in detail by van den 
Bosch & Swaters (2001), who found them to be in good 
agreement with ACDM predictions. A comparison with the 
AMDs of dark matter haloes obtained by BOO reveals the 
following: 

• The ratios of disk mass to total virial mass, /disk, in- 
ferred from mass models fitted to the observed rotation 
curves, are much smaller than the Universal baryon frac- 
tion /bar for currently popular cosmologies. This indicates 
that large fractions of baryonic mass have either been pre- 
vented from cooling, or have been removed from the disk or 
halo by either feedback or stripping mechanisms. 

• The distribution of Adisk is in good agreement with 
the distribution of halo spin parameters. This suggests that 
disks have the same total specific angular momentum as dark 
matter haloes. 

• The disk mass fractions /disk are strongly correlated 
with the disk's spin parameter Adisk- 

• The distribution of f = jtot/jmax of disks is narrower 
than p(() of dark matter haloes, and with a mean that is sig- 
nificantly higher. This implies that the AMDs of haloes have 
more extended tails to high specific angular momentum. 

• The normalized angular momentum distributions of 
(low mass) disk galaxies are clearly distinct from those of 
dark matter haloes. The latter have AMDs that extent to 
higher specific angular momentum, and with much more 
mass at low angular momentum. 
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Despite uncertainties in the disk's density distribution, re- 
lated to the unknown mass-to-light ratio and the unknown 
amount of molecular gas, these results are robust, and con- 
firm recent findings of Navarro & Steinmetz (2000) that ap- 
parently disks form out of a small fraction of the available 
baryonic mass, but yet manage to tap most of the available 
specific angular momentum. 

Understanding these findings within the standard 
framework of disk formation is challenging. Neither the main 
mechanisms that lead to small values of /disk, i.e., feedback, 
stripping or cooling, nor viscous processes that redistribute 
specific angular momentum, seem able to provide a mean- 
ingful explanation. Somehow these results seem to imply 
that the baryonic mass components out of which disk galax- 
ies form through cooling have angular momentum distribu- 
tions that are clearly distinct from those found by BOO. This 
implies either that some mechanism decoupled the baryons 
from the dark matter during the early stages of the forma- 
tion of galaxies (during which they acquire angular momen- 
tum from cosmological torques), or that the distributions 
found by BOO are poor descriptions of the actual AMDs of 
dark matter haloes (perhaps reflecting a different nature of 
the dark matter). It is clear that without a proper under- 
standing of (the origin of) the angular momentum distribu- 
tion of the baryonic mass component of protogalaxies our 
picture of the formation of disk galaxies is incomplete. 
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